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<p ■ Abstract 



We have studied the superfluid density and the specific heat of the x — y 
model on lattices L x L x H with L 3> H (i.e. on lattices representing a 

a . 

film geometry) using the Cluster Monte Carlo method. In the //-direction 
we applied staggered boundary conditions so that the order parameter on the 

o. 

top and bottom layers is zero, whereas periodic boundary conditions were 



applied in the L-directions. We find that the system exhibits a Kosterlitz- 
Thouless phase transition at the -ff-dependent temperature T 2 below the 



{SJ . critical temperature T\ of the bulk system. However, right at the critical 

o, 

temperature the ratio of the areal superfluid density to the critical tempera- 
ture is /^-dependent in the range of film thicknesses considered here. We do 
not find satisfactory finite-size scaling of the superfluid density with respect to 
H for the sizes of H studied. However, our numerical results can be collapsed 



onto a single curve by introducing an effective thickness Ft e ff = H + D (where 



D is a constant) into the corresponding scaling relations. We argue that the 
effective thickness depends on the type of boundary conditions. Scaling of the 
specific heat does not require an effective thickness (within error bars) and 
we find good agreement between the scaling function f% calculated from our 
Monte Carlo results, f\ calculated by renormalization group methods, and 



the experimentally determined function f\. 
64.60.Fr, 67.40.-w, 67.40.Kh 
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I. INTRODUCTION 

The theory of second order phase transitions is based on the assumption that at tempera- 
tures close to the critical temperature T A there is only one dominating length scale associated 
with the critical behavior of the system, the correlation length. Since the correlation length 
diverges as the critical temperature is approached the microscopic details of the system are 
irrelevant for its critical behavior. This intuitive picture has its foundation in the renor- 
malization group treatment of second order phase transitions. Within the renormalization 
group treatment it becomes evident that the critical behavior can be divided into different 
universality classes. Each universality class is characterized by a set of critical exponents 
which describe the singular behavior of physical quantities in terms of the reduced temper- 
ature t = T/T\ — 1, e.g. for a three-dimensional bulk system the correlation length £(t) 
diverges close to T A as £(£) = ^o^l - ". 

If the system is confined in a finite geometry (e.g. a cubic or film geometry) the singular- 
ities in the physical quantities are smoothed out or a crossover to lower-dimensional critical 
behavior takes place. Finite-size scaling theory jl| is thought to describe well the behavior 
of the system at temperatures close to T\. The intuitive idea behind the finite-size scaling 
theory is that finite-size effects can be observed when the bulk correlation length becomes of 
the order of the system size (for a film geometry this is the film thickness H). For a physical 
quantity O this statement can be expressed as follows f2|: 

° {t ' H) -/f^#^V (D 



0(t,H = oo) \£(t,H = oo) 

/ is a universal function depending only on the geometry and the boundary conditions 
applied. 

Though earlier experiments on superfluid helium films of finite thickness || seemed to 
confirm the validity of the approach outlined above, in a recent experiment Rhee, Gasparini, 
and Bishop [|],[5] showed that their data for the superfluid density of thick helium films do 
not satisfy Eq. ([[]) when the expected value v = 0.67 is used. (For a comprehensive review 



of experiments on 4 He to test the finite-size scaling theory cf. Ref. ||.) As an attempt 
to understand these discrepancies between theory and experiment renormalization group 
calculations for the standard Landau- Ginzburg free energy functional in different geometries 
with Dirichlet boundary conditions (vanishing order parameter at the boundary) have been 
carried out [pHH|] . New specific heat measurements [jTIJ and also a reanalysis ||T3[ of the old 



specific heat data |L4| show good agreement between the renormalization group calculations 
reported in |7HlO| and those data. These calculations demonstrated the important role 
played by the boundary conditions. In particular, periodic boundary conditions were shown 
to be inadequate compared to Dirichlet boundary conditions to describe the experimental 
specific heat data. The renormalization group calculations have determined the specific heat 
for that range of the scaling variable where the surface contribution to the specific heat is 
dominant |]8|- |T0| , p!5| (cf. also [|16j). Such field theoretical calculations are not available for 
the case of the superfluid density and the lack of scaling in the case of the superfluid density 
of helium films is not understood. Furthermore, new experiments on liquid A He under 
microgravity conditions are planned [|IjJ to examine the finite-size scaling properties of the 
specific heat. In order to test the renormalization group calculations and because of the 
reasons above, numerical investigations of the finite-size scaling properties of the superfluid 
density fl8| and the specific heat []T9| , ^0| of thin helium films have been carried out. In Refs. 



|18| . [r9|j , we used the x — y model with periodic boundary conditions in the direction of the 
film thickness H to compute the superfluid density and the specific heat of thin helium films. 
We demonstrated scaling with respect to the film thickness using the expected values for 
the critical exponents of the superfluid density and the specific heat, thus confirming the 
validity of the finite-size scaling theory. However, the obtained universal function for the 



specific heat does not match the experimentally determined universal function of Ref. [|T2 
indicating that periodic boundary conditions are only a poor approximation of the correct 
physical boundary conditions as was already demonstrated in Ref. |T(J . Later we employed 
staggered-spin boundary conditions in the top and bottom layers of the film which improves 
the agreement between the numerically computed scaling function and the experimentally 



determined scaling function of the specific heat p0| . 



Another example where the boundary conditions play a role in the scaling behavior comes 
from Ref. |21|] where the Villain model, which also belongs to the x — y universality class, 
was studied in a film geometry with open boundary conditions in the direction of the film 
thickness. The authors of Ref. |HJ extracted the thickness dependent critical temperature 
from the temperature dependence of the correlation length in the disordered phase and found 
for the critical exponent v the value v = 0.71(1) which is different from its value of 0.6705 
known from experiments on liquid Helium |22| . 

In this paper we intend to study the effect of staggered-spin boundary conditions 
(Dirichlet-like boundary conditions, i.e. vanishing order parameter on the film boundaries) 
on the finite-size scaling behavior of the superfluid density and the specific heat of A H e in a 
film geometry in detail. Dirichlet-like boundary conditions are believed to approximate the 
physical boundary conditions more closely [|I0| , [23|1 . Throughout our numerical calculations 
we are going to describe superfluid 4 H e near the A-critical point by another form of the 
standard Landau- Ginzburg free energy functional: the x — y model (cf. e.g. Ref. |24}| ). In 
the pseudospin notation the x — y model takes the following form: 

H = -Jj2?i-?j> (2) 

(iJ) 
where the summation is over all nearest neighbors, the two-component vector s = 
(cos 6, sin 6), and J > sets the energy scale. The angle 9 corresponds to the phase of 
the expectation value of the helium atom creation operator which is defined in a volume 
whose linear extensions are much larger than the interparticle spacing and much smaller 
than the correlation length. 

In this article we study the superfluid density which corresponds to the helicity modulus 
in the pseudospin notation and the specific heat of the x — y model in a film geometry, i.e. 
on L 2 x H lattices with L ^> H. The top and bottom layers are coupled to a static staggered 
spin configuration, playing the role of the "substrate" layers, so that the magnetization in 
these layers is exactly zero. The crucial difference between these boundary conditions and 



periodic boundary conditions is that the superfluid density develops a profile in the in- 
direction, whereas it is completely homogeneous for periodic boundary conditions (cf. also 
the magnetic profile for the Ising model in a film geometry with different boundary conditions 
in Ref. ||25| ). We applied periodic boundary conditions in the L-directions because we intend 
to take the limit L — > oo. In the temperature range where the model behaves effectively two- 
dimensionally we used the Kosterlitz-Thouless-Nelson renormalization group equations to 
compute the values for the helicity modulus in the L — > oo limit. This way we eliminated the 
L-dependence of our data for the helicity modulus and were able to extract the Kosterlitz- 
Thouless transition temperature T^ D (H) for different films. We investigated the validity 
of finite-size scaling for the superfluid density and the specific heat of such films of infinite 
planar dimension and finite H. We shall also discuss scaling of the experimental results for 
the specific heat and superfluid density with respect to the film thickness and compare the 
universal scaling functions to those obtained from our theoretical investigation. 

The rest of this paper is organized as follows. In the next section we introduce the 
physical observables defined for the x — y model and the numerical method we applied 
to carry out the calculations. In section [TTT], we discuss the finite-size scaling theory and 
boundary effects. In section [TV] we discuss our results and the last section briefly summarizes 
the work described in this paper. 

II. PHYSICAL OBSERVABLES AND MONTE CARLO METHOD 



For the x — y model on a cubic lattice the helicity modulus is defined as follows [26,27] : 



-£ nE si ^-^v^ }. ( 3 ) 

where V is the volume of the lattice, f3 = J/ksT, e^ is the unit vector in the corresponding 
bond direction, and e^ is the vector connecting the lattice sites i and j. In the following 
we will omit the vector index since we will always refer to the x-component of the helicity 

6 



modulus. Note that, because of isotropy, we have T x = T y . The connection between the 



helicity modulus and the superfluid density p s is established by the relation [ f28[ 

777 

Ps(T) = ( ¥ ) 2 T(T), (4) 

where m denotes the mass of the helium atom. The specific heat c is obtained from the 
energy E through 



P 2 
C= N 

where the energy is defined as: 



((E 2 ) - (Ef) , (5) 



E = £ (1 " Si • Sj) (6) 

and N is the number of spins contributing to the specific heat. 

The thermal averages denoted by the angular brackets are computed according to 

(O) = Z- 1 JUdOtOiejexpi-pE). (7) 

i 

0[6] denotes the dependence of the physical observable O on the configuration {#«}, the 
partition function Z is given by 

"n^exp(-/3£), (8) 

i 

The multi-dimensional integrals in the expressions (0) and (§) are computed by means of 
the Monte Carlo method using Wolffs 1-cluster algorithm |2^ . 



We computed the helicity modulus and the specific heat on L 2 x if lattices, where 
L = 20, 40, 60, 80, 100 and H = 4, 8, 12, 16, 20, 24. We applied periodic boundary conditions 
in the L-directions, whereas the first and the H-th layer are coupled to a boundary layer 
which consists of a staggered spin configuration, i.e. 

s(x,y) = (-iy+ys(l,l), (9) 

where x,y E [l,L] and label the integer coordinates of the lattice sites in a plane perpen- 
dicular to the if-direction. Thus, we have V = HL 2 a 3 with a denoting the lattice spacing 



and N = HI? . We carried out of the order of 20, 000 thermalization steps and of the order 
of 750, 000 measurements. The calculations were performed on a heterogeneous environ- 
ment of computers including Sun, IBM RS/6000 and DEC alpha AXP workstations and a 
Cray-YMP. 

III. FINITE-SIZE SCALING AND BOUNDARY EFFECTS 
A. The helicity modulus 



Here we shall discuss the finite-size scaling theory of the superfluid density. In Ref . |T8| . 



we studied the helicity modulus T for the x — y model in a film geometry with periodic 
boundary conditions in the if-direction and we have shown the following steps. In a certain 
temperature range around the bulk critical temperature T\ where the bulk correlation length 
£(T) becomes of the order of the film thickness H the quantity TH/T exhibits effectively two- 
dimensional behavior and a Kosterlitz-Thouless phase transition takes place at a temperature 
T^ D (H) < T\, The critical temperature T^ D (H) approaches T\ in the limit H — > oo as 



T™(H)=T X (1 + ^), (10) 

where for periodic boundary conditions we found that x c = —0.9965(9) using for the critical 



exponent v the experimental value v = 0.6705 |22j and for T\ the value T\ = 2.2017 pO 
The quantity TH/T is a function of the ratio H/£(T), i.e. 

mm =HtH v» ) . (11) 

The universal function $(x) has the properties J3TJ 

$(x > x c ) = 0, 

Hx c ) = -. (12) 

TV 

In the limit x — > x~ the function $(a;) can be written as: 



§(x) = - [1 + Ay/x^=x ), (13) 

7T v ' 



where for periodic boundary condition we found that A = 0.593(5). This form of the 
universal function reconciles the scaling expression (|TT| ) with the two-dimensional behavior 
[32[,i.e. asT ^T 2D (H) 

T(T,H)H 2 



T 



IX 



(i+H^i-^y (M) 



where |33| 



b(H) = AH 1 ' 2 ". (15) 

The results stated above confirm the theoretical expectations about scaling given by Am- 



begaokar et al. in Ref. [31 1 and agree with the experimental findings of Bishop and Reppi 



|34| and Rudnick [35]. In the case of periodic boundary conditions in the if -direction the 
validity of the finite-size scaling form (|ll]) can already be observed for films of thicknesses 
# = 6,8,10 0. 

If nonperiodic boundary conditions are introduced Privman argued that the general 
scaling form (|ll|) has to be altered into [[36 1 



T(T T g)JJ = i>(ttf ^) + u In H + d, (16) 

where u and d are constants depending on the boundary conditions. The Monte Carlo data 
for the helicity modulus obtained from a Monte Carlo simulation of the x — y model on cubes 
H x H x H where the spins in the boundary layers were all parallel (pinned-surface-spin 
boundary conditions), were found to be consistent with the presence of the logarithmic term 
in the scaling form ([TB] ) [37|. 



Let us now investigate the consequences of the logarithmic term in ([TBI) . Again we 
have to reconcile the two-dimensional behavior flT3|) and the general scaling form (|T6|) for 
temperatures close to the critical temperature T 2D (H) for thick enough films. Introducing 
the expression ( |T0"P for the if -dependent critical temperature into Eq. ( ]Tj] ) and keeping only 
terms up to H l l u under the square root leads to 

T(T,H)H 2 



T 71 



- + b{H)H- l/2u Jx c -tH l /v. (17) 



where we have absorbed the factor of 2/n in the definition of b(H). With the assumption 
(|T5|) we may make the identification 

$0) = A^x c - x. (18) 

In order to account for the logarithmic term in ([16]) we have to abandon the universal jump 
at T™(H) 

T(T™(H),H)H 2 

T? D (H) tt' l ' 

instead we have to assume 

i.e. the jump becomes if-dependent. (The numerical values for x c and A will be differ- 
ent from the values given above as they depend on the boundary conditions.) Eq. (|2"0D 
means that for nonperiodic boundary conditions in the if-direction expression ( |Tj| ) has to 
be generalized to 



^^ = 9(H) + b{H)^-^- y (21) 



where 



b(H) = AH l ' 2u , (22) 

g(H) = d + uj\xiH. (23) 

B. The specific heat 

For the finite-size scaling of the specific heat one can use similar scaling expressions to 
Eq. ([!]) which were examined in detail in Ref . |Tj| . The finite-size scaling expression for the 



specific heat c can also be written in an equivalent way as |7|||: 

(c(t, H) - c(to, oo)) E~ ajv = h{tH l l»). (24) 
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The function f\{x) is universal and v = 0.6705. At the reduced temperature to the correla- 
tion length is equal to the film thickness H, i.e. t = (^/H) 1/v with $ = 0.498 [|38j. This 
scaling form has been used to analyze the experimental data and, thus, we shall also discuss 
the scaling of the specific heat using this form in order to compare to published experimental 
results for the universal function / x . We have 

c(t ,oo) = c(0,oo) + c+UM , (25) 

where we have found that c(0, oo) = 30, cf = —30 |L9] and a/v = —0.0172 because of the 



hyperscaling relation a = 2 — 3za In Ref. [0| we demonstrated that our numerical results for 
the specific heat of the x — y model on a film geometry with periodic boundary conditions 
follow the finite-size scaling form ([M|) for the thicknesses as small as H = 6, 8, 10. 

IV. RESULTS 

Here we shall present our results for the superfluid density and the specific heat and our 
analysis for the case of staggered- spin boundary conditions as defined in ||. 

A. The helicity modulus 

In this section we would like to determine the values of the ratio T(T,H)H/T in the 
limit L — > oo and find estimates for the critical temperatures T^ D (H) and the parameters 
g(H) and b(H) (cf. Eq. (|2T|)). In order to do this we follow closely the procedure described 
in Refs. |8],|39[. 

Fig. [I] displays the Monte Carlo data for the helicity modulus in units of the lattice spacing 
a and the energy scale J for the film of fixed thickness H = 4. This figure demonstrates that 
staggered boundary conditions for the top and bottom layers of the film strongly suppress the 
values of the helicity modulus with respect to the case of periodic boundary conditions. As 
a consequence films with staggered boundary conditions have a smaller critical temperature 
than films with periodic boundary conditions. 
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At the temperature T = 2.1331 we computed the helicity modulus on a 60 x 60 x 20 
lattice for each layer separately and plotted the result T L (z)/J in Fig. 0, where z enumerates 
the layers. The layered helicity modulus is symmetric with respect to the middle layer where 
it reaches its maximum and decreases when the boundaries are approached. Although the 
helicity modulus T(T, H, L)/J is not the average of the quantity T L (z)/J over all layers (this 
is due to the second nonlinear term in expression @), the curve in Fig. |^ is an approximation 
to the profile the superfluid density develops in thin films. 

Let us turn now to the computation of the values for the ratio K = T/(TH) in the 
L — > oo limit. For a fixed thickness H and at temperatures T below but sufficiently close to 
the critical temperature T^ D (H) the system behaves like a two-dimensional system ||3T| , |T8 



In this regime we demonstrated [18 1 that the dimensionless ratio K obeys the Kosterlitz- 



Thouless-Nelson renormalization group equations [j32^,[4(i|,[41 



dK(T,l)_ A , 2 



dl 
dy(T,l) 



4vrV(T,0, (26) 

(2-nK- 1 (T,l))y(T,l). (27) 



dl 

In y is the chemical potential to create a single vortex, e denotes the size of the core radius 
of a vortex. These equations contain the universal jump K(T^ D (H),H) = n/2. In order 
to adjust the above equations to the possibility of an if-dependent jump of the ratio K at 
T^ D (H) we generalize equations (|26"D and (p7f) to: 



dK(T,l) _ 2 



(y 2 (T } l) } (28) 



dl 

2(l-eiT- 1 (T,Z)MT,0, (29) 



dy(T,l) _ nH ^_ 



dl 

where ( and e are if -dependent constants. After eliminating the variable y from the coupled 
system of differential equations we obtain: 

^^ = i(K(T,l) -e\nK(T,l) - C), (30) 

where C is a constant which satisfies the condition C > e(l — lne). This condition allows 
for the existence of roots of the right hand side of Eq. (|30|) . If we identify the scale / with 
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In L up to a constant we can use Eq. (|30|) to extrapolate the computed values K(T, H, L) 
obtained on lattices of finite planar dimension L to the L = oo limit. Namely, at L = oo 
the left hand side of Eq. fl50|) vanishes [[^.[1(].[11]] and K(T, H, oo) < e is given by the root 



of the right hand side of Eq. (|30|). The parameters e and C are found by fitting the Monte 
Carlo data for K(T,H,L) at a fixed H to the numerical solution to Eq. ([ft]). Table | 
contains our fitting results for the fitting parameters e and C and the values for K(T, H, oo) 
for the thicknesses H e [4, 20] and Fig. [| shows a typical fit. We were not able to explore 
the two-dimensional region for the film with H = 24 because the temperature range where 
the film behaves two-dimensional becomes very narrow and our computer resources did not 
allow an accumulation of data accurate enough to resolve this region. 

In Fig. | we plot X (T, H)H/T versus tH l / v for the thicknesses H = 12, 16, 20, 24 to check 
the validity of the scaling form ([TT|) where v = 0.6705. The data for the helicity modulus 
used in Fig. | have completely lost their L-dependence. We do not obtain a universal 
scaling curve, thus scaling according to the expression flTTD is not valid for the films with 



thicknesses up to if = 24. Therefore we will try to employ the scaling form ( |T6| ) which 
requires the knowledge of g{H). Since K~ 1 (T,H,oo) satisfies Eq. ( |2"T| ) for a fixed H and 
temperatures close enough to the critical temperature T^ D (H) we can fit the obtained results 
for K~ 1 {T 1 H ) oo) to Eq. (|21~1) and find an estimate for T^ D (H) and the parameters b(H) 



and g(H). In Table || we present our fitting results and Fig. |5| shows the fit to the data for 
T(T, H)H/T at H = 4. It is interesting to note that the if- dependence of the parameter g 
can be described by the formula 

g{H) = 0.338(19) In H -0.238(35) (31) 

for H G [4,20]. This is consistent with Privman's prediction [p5|j . In Fig. |6] we plot 
T(T,H)H/T - g(H) versus tH l l u for H = 8,12,16,20 where g(H) is given in Table 0, 
the bulk critical temperature T\ = 2.2017. Also the scaling form ( |T6| ) does not collapse our 
data points onto one universal curve. This situation is the same as the one Rhee, Gasparini 
and Bishop encountered when they tried to verify finite-size scaling for their data of the 
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superfluid density ||||. Their data of the superfluid density did not fall onto one universal 
curve when the scaling form ( |iT|) was employed, neither did the inclusion of a logarithmic 
term as in Eq. ( jIB[) help to achieve data collapse ||. 

Of course, one reason for the failure of scaling of our data of the helicity modulus ac- 
cording to the expressions (|TTJ) or flT5| ) could be that our thicknesses are still too small. 
On the other hand Rhee et al. use films of macroscopic sizes and do not confirm scaling. 
Furthermore, for films with periodic boundary conditions scaling of the helicity modulus 
occurs already for thicknesses as small as H = 6 |18 . 



Let us therefore pursue another line of thought |42| which we borrow from the mean field 
treatment of thin ferromagnetic films |43| . The reduced critical temperature of a ferromag- 
netic film t c (H) {t c (H) = 1 — T 2D (H)/T C where T c is the 3D bulk critical temperature of 



the ferromagnet) can be obtained from the following set of equations |43 |: 



H . . 

wtanu = — , (32) 
2A 

t c {H) = i^) 2 u\ (33) 

The lattice spacing is denoted by a and A is the extrapolation length A (cf. also Ref. [fH}]). 

Let us compute t c (H) in the limit H/{2\) 3> 1. From Fig. 2 of Ref. [^3j it is clear that 
u -> vr/2 in the limit H/(2\) -> oo. Thus, writing u = n/2 - e with < e < 1 Eq.(H) 
turns into 

2u\ 7Y 

tan e w e = u. (34) 



H 2 



Solving for u yields 



2(1 + 1) 

Inserting this into Eq.(R3) we obtain finally 



</ , - ... . (35) 



a 2 n 2 



^ - w^w (36) 

Thus, within the mean field treatment the reduced critical temperature scales with the 
correct critical exponent v = 0.5 but with an effective thickness H + 2A. It is interesting 
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that 2A appears in Eg. (p6|) . This means that we have to add twice the extrapolation length 
(for each side of the film) to the actual thickness. Furthermore, if the magnetization is 
suppressed close to the boundaries (A > 0) the critical temperature T c 2D (if ) is smaller than 
the 3D bulk critical temperature and we have to add 2A to the actual thickness. 

The lack of scaling of our data of the helicity modulus with the expected critical exponent 
v = 0.6705 indicates that the critical temperatures T 2D (H) do not satisfy Eq. fllPl) . Instead, 
due to the profile of the superfluid density we may expect an effective film thickness H e ff 
which enters the scaling expressions (|10|) and (|TTD . In close analogy to the mean field 
treatment of ferromagnetic films discussed in the paragraph above, we assume that H e ff = 
H + D where D is a constant. Indeed, for the film thicknesses H = 12, 16, 20 we obtain 
x c = -3.81(14) and D = 5.79(50) with v = 0.6705. In Fig. we plot T{T,H)H eff /T 
as a function of tH e lK for films with H = 12, 16, 20, 24 where v = 0.6705. The data for 
the helicity modulus collapse reasonably well onto a single curve. We can understand the 
increment Dasa scaling correction which renders the scaling relations (|T0| ) and (JTTJ) valid 



even for very thin films. For large thicknesses H the increment D can be neglected and 
we recover the conventional scaling forms. Of course, it is possible to invent scaling forms 
different from the structure ( |3"6"D which yield the conventional scaling expressions in the limit 
H — > oo. For example, we have obtained similarly good fitting results using the expression 
t c (H) = a\H~ x l v + a 2 H~ 2 which is also motivated by mean field theory (cf. eg. Ref. [fl8|l ). 
However, since we find Eq. (|36|) physically appealing we continue to describe the effects of 
boundary conditions by an effective thickness as we have done in this paragraph. 

In order to test further the assumption that the boundaries introduce an effective thick- 



ness into the scaling expression (|Tl~D we try to describe the thickness dependence of the 
Kosterlitz-Thouless transition temperature of the Villain model with open boundary con- 
ditions (interactions of the top and bottom layer only with the interior film layers) |^TJ by 
Eq. (|T0"D, where H is replaced by the effective thickness H e ff = H + Dy. Indeed, taking 
v = 0.6705 we find 
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( T™(H) Y 1.384(9) 

{ T x ) "if + 1.05(2)' ld7j 

thus ZV = 1.05(2) and x c = —1.62(2). The function (J37|) is the solid line in Fig. |j. Again 
for this case, the increment D is a correction which makes the scaling relations ( |10D and ((TTJ) 
valid even for very thin films. The result (EH1) means that the film thicknesses considered in 



Ref. [21] were still too small to extract the expected value of the critical exponent v from 
the if -dependence of the critical temperature ( |10D without the help of an effective thickness 
H + D v . 

In Fig. ^| we achieve approximate collapse of the experimental data of Rhee et al. [|],|5| for 
the superfluid density p s for films of various thickness d (d is in /im) by plotting p s (t, d)d e ff/p 
versus tdj^ with v = 0.6705 and <i e // — d+ 0.145. This value of the effective thickness 
was found by examining the reduced temperatures tf s (H) where finite-size effects set in. 
According to finite-size scaling theory tf s has to fulfill the relation tf s ex d~ l l p ', thus in our 
case tf s oc d e ji v '. The data corresponding to the film with d = 3.9pm deviate from the 
universal curve and we attribute this to the anomalous behavior of these data. Namely, in 
general \tf s (di)\ > |£/ s (<i2)| if d\ < d 2 , but this is not the case for d\ = 2.8pm and d 2 = 3.9pm 
(cf. Refs. §||). 

Let us compare the increments over the film thickness for the three cases of film geometry 
considered above. For the Villain model with open boundary conditions we obtained Dy = 
1.05(2) while for the x—y model with staggered boundary conditions we obtained D = 5.8(5). 
All increments are expressed in lattice spacing units. The difference in these values of the 
increments reflect how severe the effect of the boundary conditions is. Open boundary 
conditions are less demanding on the order parameter at the boundary compared to staggered 
boundary conditions used in our simulations. The value of the increment in the case of 4 He 



on silicon is large, d e jf—d = 491.5 in lattice spacing units a (a = 2.95A [T9J). This might cast 
some doubts on our proposed scaling form for the superfluid data of Rhee et al. It is possible, 
however, to imagine that on the surface of these films vortices are pinned by impurities or 
other forms of disorder, which make the order parameter vanish at the boundary and which 
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introduce an effective length scale of such a magnitude. 

So far we have seen that different boundary conditions create different effective thick- 
nesses. The influence of the boundaries vanishes for thick enough films and only in a certain 
small range of film thickness the influence of the boundary conditions has to be taken into 
account. According to our findings the scaling form ([T^) for the helicity modulus in the 
limit T —* T^ D (H) takes the following form now: 

T(T, H)H eff _ ( f 1/v 



-g\l + A^jx c - tH%)\ , (38) 

where A, x c and g are constants depending on the boundary conditions. Especially for 
tH e ff = x c we should have 

T(T^H)H efJ = _ (x) 

1 c 

and this value of g is not necessarily equal to 2/tt. In Table [TH] we give the values for g found 
by our fitting procedure. We still have a slight if-dependence in g but for the thicknesses H > 
16 the value for g seems to saturate at g(H — > oo) ~ 0.97. Since we can neglect the effective 
thickness for very large film thicknesses, this means that films with staggered boundary 
conditions in the if -direction of the film exhibit a jump in the quantity T(T^ D ,H)H/T^ D 
that is different from 2/-7T which was found for films with periodic boundary conditions 



34|40lJlgl ■ Therefore, assuming that our extrapolation to large film thicknesses from small 
size films using the idea of the effective thickness is valid, we have to conclude that the jump 
T(T 6 2D , H)H/T^ D depends on the boundary conditions. In principle there is nothing wrong 
with this conclusion because the universal functions (and the jump is a particular feature 
of a particular universal function) depend very importantly (especially near the critical 
temperature) on the boundary conditions. The scaling function should not be confused 
with the critical exponents which are independent of the geometry and boundary conditions. 
The scaling functions for given universality, given geometry and boundary conditions are 
universal. This leads us to the conclusion that this jump in the experimental findings 



should depend on the substrates (cf. also Ref. p5|). This influence of the substrate must be 
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mediated by the vortices whose generation is enhanced close to the boundaries due to the 
effect of the boundary (for a more detailed discussion cf. section [V]). In the experiments 
the value of 2 ji\ was found p4| , |35| , however, the thicknesses of these films are much smaller 



than the above length scale D found to fit the data of Rhee et al. We believe that the 
vortex density was almost homogeneous throughout the films used for these measurements. 
This situation corresponds to the x — y model in a film geometry with periodic boundary 
conditions along the film direction where the vortex density is the same everywhere. 

B. The specific heat 

In this section we would like to investigate the finite-size behavior of the specific heat 
c(T, H). Since we do not possess an easily handable procedure to take the L — > oo limit for 
the values of the specific heat c(T, H, L) computed on finite lattices LxLxH we approximate 
films with infinite planar dimension by 100 x 100 x H lattices. This seems justified because 
the specific heat appears independent of L for L > 60 (cf. Fig. [10]). Furthermore, we do 
not expect the maximum of the specific heat to grow dramatically with increasing values 
of L because for temperatures in the range T^ D (H) <T<T\ the behavior of the specific 
heat can be described by the Kosterlitz-Thouless theory which leads to a finite value of this 



maximum. In order to illustrate this argument we show in Fig. [LT] the size dependence of 
the specific heat c(T, L) computed on pure two-dimensional lattices L x L with periodic 
boundary conditions. The L-dependence of the specific heat can be neglected for values of 
L>80. 



In Fig. 12 we compare the specific heat of films with H = 16, 20, 24 to the bulk specific 



heat c(t,H = oo) taken from Ref. fl~9 |. The specific heat values for films of thickness 



H = 12, 16, 20 lie above the bulk curve for t < 0. Such a behavior is also found in experiments 
on helium films about 30A thick |J. This crossing effect is due to the large shift of the 
temperature T m where the specific heat for a certain film thickness takes its maximum down 
to temperatures below the bulk critical temperature T\. For thicker and thicker films the 



maximum temperature T m approaches T\ , thus the confined specific heat data will fall 
below the bulk curve which is expected from the field theoretical calculations P|- ]ID , ^5| . The 



film thicknesses used in our Monte Carlo calculation range from H = 35A to H = 70. 8A 
(the lattice spacing a = 2.95A |0|) which is comparable to the experimental film thicknesses 
(30A) where this crossing effect can be observed. Fig. [IB] shows the specific heat of the film 
with H = 24 alone, indicating that the effect of crossing the bulk curve indeed vanishes for 
thicker and thicker films. Unfortunately it is beyond our means to carry out the necessary 
analysis for thicker films than were treated here. It is interesting to note that we find the 
same qualitative behavior of the specific heat in the case of the x — y model in a cylindrical 
geometry 



In Fig. [14| we plot the scaling function fi(x) given by expression (|24]). According to the 
previous discussion this scaling function can only be an approximation to the correct one 
which one would need to compute from films with H > 24 and L » H . We find that 
our data for the specific heat for films of various thicknesses collapse approximately onto a 
single curve. It seems that the specific heat is rather insensitive to the boundary effect of 
introducing an effective thickness and a very high accuracy in the computation of the specific 
heat is needed to detect the presence of the effective thickness in the scaling function fi(x). 
For example, one could determine the temperatures T m (H) where the specific heat reaches 
its maximum and examine the validity of Eq. (|10"D for T m (H), because the if -dependence 
of T m (H) is also given by expression fllTf ) (with a different value for x c than for the critical 
temperatures T^ D (H)). This can be done more easily in experiments because in Monte 
Carlo calculations an extrapolation procedure for the values of the specific heat at finite 
planar dimensions to the values at infinite planar dimensions is needed (and which is not 
available at present) whereas the films used in experiments represent films with infinite 
planar dimensions. 

We can directly compare our function fi{tH l ^ u ) to the experimentally determined scaling 
function f\{x) given in Refs. [|I2[] by expressing all lattice units in physical units using the 



conversion formula: 
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-, / m V m k B ( a\ a,v , , Nl .___ Joule 

/iWlp^ = "^T" (Xj /l( X )l^ C e = 15 - 02 o K mQle AOOIjottfce , ( 4 0) 

where V m is the molar volume of 4 if e at saturated vapor pressure at T\, the lattice spacing 
a = 2.95A |TP| and the film thickness is measured in A. In Fig. [T5| we compare the functions 
fi(x) obtained from Monte Carlo calculations of the specific heat of films with periodic 
boundary conditions and staggered boundary conditions in the direction of the film thickness 
to the experimentally determined function f\{x) and to the function fi(x) obtained from 
field theoretical calculations j|,|9| . This figure clearly shows the influence of the boundary 
conditions on the shape of the universal function as was already demonstrated by the field 
theoretical calculations reported in Ref. [IIJ. In Fig. [15| we see that the scaling function 



fi(x) for films with staggered boundary conditions crosses the scaling function fi(x) for 
films with periodic boundary conditions (cf. also p0| ) in the range — 14 < tH 1 ^ < —8 (cf. 
Fig. |15|) with H measured in A. This crossing is due to the relative smallness of our film 
thicknesses (see the discussion above) and does not occur if we had used much thicker films 
in our Monte Carlo calculations to deduce the scaling function fi(x) pBfl . We expect our 
function fi(x) to be slightly modified in the range — 14 < tH 1 ^ < —8 (cf. Fig.|l|) when 
it is computed from much thicker films. We believe that the wings of our curve, however, 
will remain unchanged. Unfortunately, it is beyond our computational means to repeat the 
calculations for thicknesses larger than 24. 

V. DISCUSSION OF THE RESULTS 

Our findings suggest that it is possible to introduce an effective thickness H e ff = H + D 
into the scaling expressions for the superfluid density (|Tl"| ) and achieve scaling (i.e. data 
collapse) even for rather thin films. The increment D over the film thickness H can be 
understood as an effective correction to scaling. The appearance of the effective thickness in 
our scaling-function can be understood as follows. The superfluid density has a finite value 
in the first layers next to the boundary layers of the film, and its rise from this finite value 
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to its bulk value p™(T, H) (for T < T^ D (H)) inside the film can be divided into two regions. 
Let us assume that T is very close to T^ D (H) where the correlation length is very large. 
There is a rather narrow region of thickness D\ of the film which is in contact with the 
boundary wall where the superfluid density rises very fast to attain some value pf 1 < p™. 
Then it rises with a much slower rate over a length scale of the order of the correlation 
length to reach its value of p™. The reason for the initial fast rise are the correlations over 
length scales much smaller than the correlation length. One might think that one then 
has to exclude the region of the film where the superfluid density rises very sharply and 
this leads to a negative value for D. However, this initial rise of p s occurs very fast over 
a length scale Di which is smaller than a length D 2 which would have been required in 
order for the superfluid density to reach the same value if this rise would have occurred over 
larger distances over which the spin-spin correlations are governed by the correlation length 
which controls the long distance behavior of the correlation function. This implies that the 
required increment to the thickness is D = D 2 — -Di, which is positive. 

Only for films with thicknesses H which fulfill H 3> D scaling with H can be observed. 
For periodic boundary conditions we have D = [TE\, open boundary conditions seem to 



yield D = 1.05 and for staggered boundary conditions we obtain D = 5.79. Due to their 
structure staggered boundary conditions support vortex formation close to the boundaries, 
thus the superfluid density decreases from its maximum in the middle of the film towards 
the boundaries. This effect is less pronounced for open and absent for periodic boundary 
conditions. Our results imply that the more the superfluid density is suppressed near the 
boundaries the larger is the value of D. The data for the superfluid density of Rhee et al. 
||] which correspond to 4 He on Si require a large value of D = 0.145pm. Thus, Si should 
suppress the superfluid density dramatically close to the boundary. Since D is so large, only 
films with H > 3pm should allow scaling with H. It would be interesting to investigate the 
scaling behavior of the superfluid density of 4 He on different substrates (which represent 
different types of boundary conditions) in a wide range of film thicknesses to check our 
hypothesis. 
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A consequence of scaling our data of the helicity modulus (or superfluid density) using 
an effective thickness is that the jump in the quantity T(T^ D ,H)H/T^ D depends on the 
boundary conditions in the top and bottom layers of the film, i.e. on the substrates in real 
helium experiments. The value of the jump is only 2/n (in lattice units) for film thicknesses 
small compared to D as is the case in the experiments reported in Refs. [|51|,ft5| . The same 
value of the jump was found is the case for x — y films with periodic boundary conditions 
|18| . Thus, experiments could be also used to determine the jump in the areal superfluid 



density at the Kosterlitz-Thouless temperature and determine its substrate dependence. 

In this work we represented the substrate by a staggered spin configuration coupled to 
the top and bottom layers of the film to simulate Dirichlet-like boundary conditions (zero 
order parameter in the boundary) in the substrate. For the staggered spin configuration 
the local magnetization is exactly zero on a plaquette, i.e. on a domain of the size a x a 
where a denotes the lattice spacing. Dirichlet-like boundary conditions in the substrate are, 
however, also realized by the following spin configuration: The spins are parallel in a square 
domain of linear dimension R, but any two spins representing two adjacent domains are 
antiparallel to each other. Thus, an additional length scale associated with a finite value 
of the local magnetization over the length scale R is introduced and influences the scaling 
behavior of the helicity modulus and the specific heat. Since disorder in the boundaries 
supports vortex formation close to the boundaries, vortices should play an active role in 
creating the boundary effect described above. 

VI. SUMMARY 

We have investigated the finite-size scaling properties of the specific heat c and the 
helicity modulus T of the x — y model in a film geometry, i.e. on L x L x H lattices 
with L ^> H where staggered and periodic boundary conditions where applied in the in- 
direction and the L-directions of the film, respectively. We found that a Kosterlitz-Thouless 
phase transition takes place at the if-dependent critical temperatures T c 2D (if), however, 
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for the films used in our calculations the jump T(T C 2 (H), H)H/T^ (H) appears to be 



//-dependent. Furthermore, scaling of the helicity modulus according to Eq. ( |TT| ) is not 



valid for our film thicknesses, neither is scaling according to Eq. (|T6| ) which was derived 



following a suggestion of Privman |56| . Introducing an effective thickness H e ff > H into 
the scaling expression (JTTj) we are able to collapse our data as well as the data of Rhee 
et al. Q reasonably well onto one universal curve. Our results suggest that the boundary 
effect of creating an effective thickness H e ff depends on the boundary conditions which can 
be realized in experiments by different substrates and is negligible for thicknesses H which 
fulfill H eff /H -1«1. We argue that the jump in the quantity T(T 2D , H)H/T? D depends 
on the boundary conditions and is 2/tt only for certain ideal boundary conditions such as 
the periodic boundary conditions. Within error bars scaling of the specific heat does not 
require an effective thickness and the scaling function / x (x) for the specific heat agrees rather 
well with the experimentally determined scaling function fi(x) and with the result of the 
renormalization group calculations reported in [§,0. However, we found that Monte Carlo 
simulations of much thicker films than we have used have to be performed to determine the 
position of the maximum of the scaling function accurately. At present this is unfortunately 
beyond our computer resources. 

VII. ACKNOWLEDGEMENTS 

We would like to thank Prof. Dohm for interesting and useful discussions. N. S. would 
like to thank the Hochsleistungsrechenzentrum Julich for the opportunity of using their 
computing facilities. This work was supported by the National Aeronautics and Space Ad- 
ministration under grant no. NAG3-1841 and by Sonderforschungsbereich 341 der Deutschen 
Forschungsgemeinschaft . 



23 



REFERENCES 

[1] M. E. Fisher and M. N. Barber, Phys. Rev. Lett. 28, 1516 (1972); M. E. Fisher, Rev. 
Mod. Phys. 46, 597 (1974); V. Privman, Finite Size Scaling and Numerical Simulation 
of Statistical systems, Singapore: World Scientific 1990. 



[2] 
[3] 
[4] 
[5] 

[6] 
[7] 



[9 
[10 

[11 
[12 

[13 
[14 

[15 
[16 
[17 



E. Brezin, J. Physique 43, 15 (1982). 

J. Maps and R. B. Hallock, Phys. Rev. Lett 47, 1533 (1981). 

I. Rhee, F. M. Gasparini, and D. J. Bishop, Phys. Rev. Lett. 63, 410 (1989). 

I. Rhee, D. J. Bishop, and F. M. Gasparini, Physica B165&166, 535 (1990). 

F. M. Gasparini and I. Rhee, Prog. Low Temp. Phys. XIII, 1 (1992). 

R. Schmolke, A. Wacker, V. Dohm, and D. Frank, Physica B165 & 166, 575 (1990). 

V. Dohm, Physica Scripta T49 46 (1993). 

P. Sutter and V. Dohm, Physica B194-196, 613 (1994). 

W. Huhn and V. Dohm, Phys. Rev. Lett. 61, 1368 (1988). 

M. Krech and S. Dietrich, Phys. Rev. A46 1886 (1992), 1922 (1992). 

J. A. Nissen, T. C. P. Chui, and J. A. Lipa, J. Low Temp. Phys. 92, 353 (1993), Physica 
B194-196, 615 (1994). 

A. Wacker and V. Dohm, Physica B194-196 611 (1994). 

T. Chen and F. M. Gasparini, Phys. Rev. Lett. 40, 331 (1978); F. M. Gasparini, T. 
Chen, and B. Bhattacharyya, Phys. Rev. 23, 5797 (1981). 

M. Pitz, Diploma thesis, RWTH Aachen, Germany (1996). 

M. Coleman and J. A. Lipa, Czech. J. Phys. 46 (Suppl. SI), 183 (1996). 

J. A. Lipa, private communications. 

24 



[18] N. Schultka and E. Manousakis, Phys. Rev. B51, 11712 (1995). 

[19] N. Schultka and E. Manousakis, Phys. Rev. B52, 7528 (1995). 

[20] N. Schultka and E. Manousakis, Phys. Rev. Lett. 75, 2710 (1995). 

[21] W. Janke and K. Nather, Phys. Rev. B48 15807 (1993). 

[22] L. S. Goldner and G. Ahlers, Phys. Rev. B45, 13129 (1992). 

[23] V. L. Ginzburg and L. P. Pitaevskii, Soviet Physics, J. E. T. P., 34, 858 (1958). V. L. 
Ginzburg, Soviet Physics, J. E. T. P. 2, 170 (1956). 

[24] H. Kleinert, Gauge Fields in Condensed Matter, Singapore: World Scientific 1989. 

[25] S. Dasgupta, D. Stauffer, and V. Dohm, Physica A213, 368 (1995). 

[26] S. Teitel and C. Jayaprakash, Phys. Rev. B27, 598 (1983). 

[27] Y.-H. Li and S. Teitel, Phys. Rev. B40, 9122 (1989). 

[28] M. E. Fisher, M. N. Barber and D. Jasnov, Phys. Rev. B16, 2032 (1977). 

[29] U. Wolff, Phys. Rev. Lett. 62, 361 (1989). 

[30] W. Janke, Phys. Lett. A148, 306 (1992). 

[31] V. Ambegaokar, B. I. Halperin, D. R. Nelson and E. D. Siggia, Phys. Rev. B21 1806 
(1980). 

[32] D. R. Nelson and J. M. Kosterlitz, Phys. Rev. Lett. 39, 1201 (1977). 

[33] R. G. Petschek, Phys. Rev. Lett. 57, 501 (1986). 

[34] D. J. Bishop and J. D. Reppy, Phys. Rev. Lett. 40, 1727 (1978). 

[35] I. Rudnick, Phys. Rev. Lett. 40, 1454 (1978). 

[36] V. Privman, J. Phys. A23, L711 (1990). 



25 



[37] K. K. Mon, Phys. Rev. B44, 2643 (1991). 

[38] A. P. Gottlob and M. Hasenbusch, Physica A201, 593 (1993). 

[39] N. Schultka and E. Manousakis, Phys. Rev. B49, 12071 (1994). 

[40] J. V. Jose, L. P. Kadanoff, S. Kirkpatrick, and D. R. Nelson, Phys. Rev. B16, 1217 
(1977) 

[41] J. M. Kosterlitz, J. Phys. C7, 1046 (1974); J. M. Kosterlitz and D. J. Thouless, J. Phys. 
C6, 1181 (1973). 

[42] N. Schultka and E. Manousakis, J. Low Temp. Phys. 105, 3 (1996). 

[43] M. I. Kaganov and A. N. Omel'yanchuk, Sov. Phys. JETP 34, 895 (1972). 

[44] K. Binder and P. C. Hohenberg, Phys. Rev. B 6, 3461 (1972). 

[45] N. Schultka and E. Manousakis, Czech. J. Phys. 46 (Suppl. SI), 451 (1996). 

[46] V. Dohm, private communication. 

[47] N. Schultka and E. Manousakis, in preparation. 

[48] J. Rudnick and G. Gaspari, Phys. Rev. B32, 7594 (1985). 



26 



FIGURES 



0.5 



0.4 



3 0.3 



X 



k 




x 40x40x4 pbc 
\ 60x60x4 pbc 
o 100x100x4 pbc 
v 20x20x4 sbc 
* 40x40x4 sbc 
□ 60x60x4 sbc 
c 100x100x4 sbc 




FIG. 1. The helicity modulus T(T, H, L) as a function of T for various lattices I? x 4 with 
staggered boundary conditions (sbc) and periodic boundary conditions (pbc) in the if -direction. 
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FIG. 2. The approximate profile T^(z) of the helicity modulus computed on a 60 x 60 x 20 
lattice at T = 2.1331, i.e. close to the critical temperature, Tj? (20) = 2.1346. 
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FIG. 3. T/(T(T,H,L)H) as a function of InL at T = 1.8182 and H = 4. The solid curve is 
the fit to the solution to (pQ). 
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FIG. 4. T(T,H)H/T as a function of tH 1 ' u for various thicknesses, v = 0.6705. 
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FIG. 5. T(T, H)H/T at L = oo and H = 4 as a function of T. The solid curve is the fit to (21) 
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FIG. 6. T(T, H)H/T - g{H) as a function of tH 1 ^. u = 0.6705 and the values for g{H) are 
taken from Table ITTL 
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FIG. 7. T(T, H)H e ff/T as a function of tHjJi for various thicknesses. -He// = H + 5.79 and 



i/ = 0.6705. 
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FIG. 8. (1 — T% (H) /T\) v for the Villain model in a film geometry with open boundary 
conditions as a function of H. The solid line is the expression (| 
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FIG. 9. Scaling of the superfluid density data of Rhee et al. [|j with the effective thickness 
d e ff = d + 0.145. v = 0.6705 and all lengths are in fim. 
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FIG. 10. The specific heat c(T, H, L) as a function of T for I? x 12 lattices. T c 2 -°(12) = 2.086, 
T x = 2.2017. 
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FIG. 11. The specific heat c(T,L) for pure two-dimensional lattices L x L with periodic 
boundary conditions. 
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FIG. 12. The specific heat c(t, H) for various films of finite thickness H compared to the bulk 
specific heat c(t, oo). 
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FIG. 13. The specific heat c(t, H) for the film of thickness H = 24 compared to the bulk 
specific heat c(t, oo). 
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FIG. 14. Scaling function /i(x) for films with staggered boundary conditions where x = tH 1 '". 
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FIG. 15. The experimentally determined function fi(x) (open circles) of Refs. [O], fi(x) for 
films with periodic boundary conditions (triangles) of Ref. |0| and staggered boundary conditions 
(filled circles), and f\(x) for films with Dirichlet boundary conditions of Ref. |g] (dashed line). 
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TABLES 



H 


T 


e 


C 


K(T,H, oo) 


x 2 


Q 


4 


1.8315 
1.8305 
1.8298 
1.8290 
1.8265 
1.8248 
1.8198 
1.8182 
1.8149 
1.8116 
1.8083 


3.2080(6) 

3.2251(4) 

3.2222(5) 

3.2125(5) 

3.0255(7) 

3.0353(8) 

3.1442(6) 

3.0890(7) 

3.5845(1) 

5.21092(2) 

10.53880(4) 


-0.51(1) 

-0.511(9) 

-0.495(9) 

-0.470(9) 

-0.262(7) 

-0.243(7) 

-0.244(6) 

-0.186(3) 

-0.478(8) 

-1.46(1) 

-4.60(4) 


2.852(87) 

2.742(50) 

2.674(43) 

2.601(37) 

2.453(29) 

2.347(23) 

2.122(12) 

2.0866(72) 

1.9916(9) 

1.9094(7) 

1.8434(8) 


0.74 

0.16 

0.17 

0.15 

0.58 

0.45 

0.02 

0.66 

0.03 

0.009 

0.06 


0.39 
0.69 
0.68 
0.70 
0.45 
0.50 
0.89 
0.52 
0.86 
0.92 
0.81 


8 


2.0167 
2.0161 
2.0155 
2.0147 
2.0141 
2.0135 
2.0127 
2.0121 
2.0115 
2.0107 
2.0101 
2.0094 


2.0420(2) 

2.0787(2) 

2.1046(2) 

2.0465(1) 

2.0371(1) 

2.0262(1) 

2.41711(9) 

2.50265(9) 

2.53788(9) 

4.78735(2) 

5.54710(2) 

3.39908(6) 


0.624(3) 

0.614(3) 

0.608(3) 

0.641(2) 

0.650(2) 

0.659(2) 

0.509(3) 

0.483(4) 

0.477(4) 

-0.383(9) 

-0.65(1) 

0.183(6) 


1.666(14) 

1.633(12) 

1.608(11) 

1.5890(84) 

1.5750(82) 

1.5610(83) 

1.5190(57) 

1.5032(53) 

1.4908(57) 

1.4739(42) 

1.4629(42) 

1.4528(45) 


1.76 
1.63 
1.58 
0.16 
0.12 
0.07 
3.12 
3.20 
3.08 
1.56 
1.48 
1.44 


0.18 
0.20 
0.21 
0.69 
0.73 
0.79 
0.08 
0.07 
0.08 
0.21 
0.22 
0.23 


12 


2.0846 
2.0842 


1.7021(3) 
1.7190(2) 


0.818(3) 
0.817(3) 


1.446(14) 
1.420(12) 


1.17 

0.77 


0.28 
0.38 
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2.0838 
2.0833 
2.0829 
2.0825 
2.0822 
2.0818 
2.0812 
2.0805 


1.7148(3) 
1.7221(1) 
1.7664(2) 
1.7455(2) 
1.7343(2) 
2.1259(1) 
2.1171(1) 
2.1161(1) 


0.822(3) 
0.822(2) 
0.817(2) 
0.826(2) 
0.831(2) 
0.733(3) 
0.745(3) 
0.754(3) 


1.406(12) 

1.3949(86) 

1.3582(64) 

1.3500(66) 

1.3466(68) 

1.3140(48) 

1.2998(47) 

1.2859(49) 


0.52 

0.25 

1.34 

1.32 

1.29 

0.004 

0.007 

0.06 


0.47 
0.78 
0.26 
0.27 
0.28 
0.95 
0.93 
0.81 


16 


2.1173 
2.1171 
2.1169 
2.1164 
2.1160 
2.1153 
2.1148 
2.1142 
2.1119 


1.4971(3) 
1.5517(3) 
1.5591(3) 
1.5589(3) 
1.7307(2) 
2.0342(4) 
1.7092(3) 
1.5761(3) 
1.6254(5) 


0.903(3) 
0.894(3) 
0.895(3) 
0.904(3) 
0.880(3) 
0.842(3) 
0.896(2) 
0.921(2) 
0.936(4) 


1.327(21) 

1.294(14) 

1.282(13) 

1.241(12) 

1.2104(61) 

1.1825(48) 

1.1816(56) 

1.1748(68) 

1.1197(78) 


1.36 
1.56 
1.53 
0.14 
1.64 
1.05 
2.20 
0.40 
0.42 


0.24 
0.21 
0.22 
0.87 
0.19 
0.35 
0.11 
0.67 
0.66 


20 


2.1336 
2.1331 
2.1327 
2.1322 
2.1317 
2.1313 
2.1308 
2.1299 


1.280(1) 
1.273(2) 
1.273(3) 
1.238(2) 
1.310(3) 
1.303(2) 
1.293(1) 
1.692(2) 


0.972(2) 
0.974(2) 
0.977(2) 
0.981(2) 
0.982(2) 
0.986(2) 
0.989(2) 
1.00(1) 


1.146(18) 

1.130(16) 

1.109(14) 

1.111(17) 

1.066(10) 

1.0535(94) 

1.0428(91) 

1.001(14) 


0.32 
0.45 
0.55 
0.13 
0.05 
0.04 
0.06 
0.03 


0.73 
0.64 
0.58 
0.88 
0.95 
0.96 
0.94 
0.86 
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TABLE I. Fitted values of the parameters e and C of the expression (30) and the extrapolated 
values K (T, H, oo) at various temperatures T and different thicknesses H. x 2 and the goodness of 
the fit Q are also given. 



H 


g(H) = r(T™(H),H)H/T? D (H) 


b(H) 


T 2 C D {H) 


x 2 


Q 


4 


0.231(13) 


2.564(63) 


1.8354(9) 


1.15 


0.33 


8 


0.47(10) 


2.94(83) 


2.0207(41) 


0.88 


0.54 


12 


0.587(44) 


3.68(57) 


2.0862(10) 


1.32 


0.24 


16 


0.715(32) 


4.07(78) 


2.1175(3) 


0.23 


0.87 


20 


0.754(53) 


4.91(91) 


2.1346(7) 


0.73 


0.60 



TABLE II. Fitted values of the parameters g(H), b(H) and TJ? (H) of the expression (Ell) for 
different thicknesses H . x 2 an d the goodness of the fit Q are also given. 



2D i 



g = T(T c w (H),H)H eff /T™(H) 



0.565(32) 

0.81(17) 

0.870(65) 

0.974(44) 

0.972(68) 



TABLE III. The jump g = T(T^ U (H), H)H eff /T^ u (H) for different thicknesses H. 
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